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A  Physics-Based  Frequency  Dispersion 
Model  of  GaN  MESFETs 

Syed  S.  Islam,  Member,  IEEE,  A.  F.  M.  Anwar,  Senior  Member,  IEEE,  and  Richard  T.  Webster,  Member,  IEEE 


Abstract — A  physics-based  model  for  GaN  MESFETs  is 
developed  to  determine  the  frequency  dispersion  of  output  resis¬ 
tance  and  transconductance  due  to  traps.  The  equivalent  circuit 
paranieters  are  obtained  by  considering  the  physical  mechanisms 
for  current  collapse  and  the  associated  trap  dynamics.  Detrapping 
time  extracted  from  drain-lag  measurements  are  1.55  and  58.42 
s  indicating  trap  levels  at  0.69  and  0.79  eV,  respectively.  The 
dispersion  frequency  is  in  the  range  of  megahertz  at  elevated 
temperature,  where  a  typical  GaN  power  device  may  operate, 
although  at  room  temperature  it  may  be  few  hertz.  For  a  1.5  x 
150  fim  GaN  MESFET  with  drain  and  gate  biases  of  10  V  and 
—IV,  respectively,  5%  decrease  in  transconductance  and  62% 
decrease  in  output  resistance  at  radio  frequencies  (RFs)  from  their 
DC  values  are  observed.  The  dispersion  characteristics  are  found 
to  be  bias  dependent  A  significant  decrease  in  transconductance  is 
observed  when  the  device  operates  in  the  region  where  detrapping 
is  significant.  As  gate  bias  approaches  toward  cutoff,  the  difference 
between  output  resistance  at  dc  and  that  at  RF  increases.  For  drain 
and  gate  biases  of  10  and  —5  V,  output  resistance  decreases  from 
60.2  kfl  at  dc  to  7.5  kfi  at  RF  for  a  1.5  /xm  X  150  GaN  MESFET. 

Index  Terms — Current  collapse,  detrapping  time,  frequency  dis¬ 
persion,  GaN  MESFET. 


I.  Introduction 

GaN-BASED  heterojunction  field  effect  transistors  (FETs) 
that  offer  a  high  bandgap  energy  of  3.4  eV  and  a  high 
breakdown  electric  field  of  2.7  x  10^  V/cm  [1]  are  highly 
attractive  for  high-power  and  high-temperature  applications  at 
microwave  frequencies.  Recently,  GaN  HEMTs  with  a  power 
output  up  to  9.8  W/mm  at  8  GHz  [2]  have  been  demonstrated. 
With  SiC  as  a  substrate,  GaN  FETs  have  been  operated  up 
to  750  °C,  as  has  been  demonstrated  by  Daumiller  et  aL  [3]. 
With  a  high  electron  saturation  velocity  of  2.7  x  10^  cm/s 
[4]  and  low  dielectric  constant  with  a  commensurate  lower 
capacitance,  GaN-based  devices  are  potential  choices  for 
high-frequency  operation.  GaN  HEMT  with  fx  =  101  GHz 
and  /max  =  155  GHz  has  also  been  demonstrated  [5]. 

GaN-based  devices  are  plagued  with  traps  that  lead  to  current 
collapse  in  the  cureent-voltage  (/-V)  characteristics  [6]-[9], 
which  results  in  severe  performance  degradation  at  high-power 
and  high-fi:equency  applications.  Zhang  et  al  [6]  reported 
current  collapse  in  epitaxially  grown  GaN  junction  FETs  at 
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high  drain  bias  and  proposed  an  estimation  of  the  length  of  the 
trapped  charge  region  along  the  channel,  at  the  channel-sub¬ 
strate  interface.  Klein  et  al  [7]  have  reported  current  collapse 
in  GaN  MESFETs  in  the  absence  of  light.  They  also  presented 
a  theoretical  model  describing  the  optical  restoration  of  the 
collapsed  drain  current  at  low  drain  bias.  Binari  et  al  [8] 
have  reported  current  collapse  at  moderately  high  drain  bias 
in  GaN  HEMTs.  Kunihiro  et  al  [9]  have  reported  recovery  of 
current  collapse  with  measurements  with  10  s  hold  time  as  the 
electrons  get  more  time  to  be  detrapped.  At  high  frequencies, 
when  the  detrapping  time  is  comparable  to  or  greater  than 
the  signal  period,  traps  cannot  respond  quickly  enough  to  the 
applied  voltage  giving  rise  to  frequency  dispersion  in  device 
transconductance  and  output  resistance  [10],  [11]. 

For  microwave  circuit  design,  the  frequency  dispersion 
of  electrical  characteristics  is  important,  as  trapping  effects 
reduce  the  transconductance  and  output  resistance  significantly 
from  their  dc  values.  Moreover,  the  large-signal  device  models 
required  to  design  analog  and  mixed-signal  and  high-power 
circuits  must  include  dc  to  radio  frequency  (RF)  dispersion  of 
the  device  characteristics  for  accurate  analysis  [10],  [12]-[14]. 
Golio  et  al  [12]  used  an  additional  voltage-dependent  current 
source  in  an  RC  network  to  account  for  the  frequency  dispersion 
of  transconductance  and  output  resistance  in  GaAs  MESFETs. 
Ytterdal  et  al  [13]  reported  an  AlGaAs-GaAs  heterostructure 
FET  model  suitable  for  mixed-mode  circuit  design  incorporating 
frequency  dispersion  of  transconductance  and  output  resistance. 
In  these  works,  the  circuit  parameters  were  extracted  from 
experimental  data.  Lee  et  al  [10]  reported  a  physics-based 
model  to  address  frequency  dispersion  in  GaAs  MESFETs  due 
to  deep-level  traps  by  obtaining  the  self-backgating  parameters 
from  experimental  data.  Larson  [14]  proposed  an  improved 
circuit  model  to  incorporate  the  bias  dependence  of  output  resis¬ 
tance  dispersion.  As  discussed,  the  determination  of  frequency 
dispersion  of  output  resistance  and  transconductance  has  been 
based  upon  extracting  circuit  parameters  from  experimental  data 
and  not  on  physical  principles  describing  trapping  effects.  In 
this  paper,  the  applied  bias,  detrapping  time,  and  temperature 
dependences  of  frequency  dispersion  of  output  resistance  and 
transconductance  are  investigated  by  formulating  the  equivalent 
circuit  parameters  of  the  intrinsic  GaN  MESFET  obtained 
mainly  from  the  underlying  physics  of  the  device. 

n.  THEORY 

The  effect  of  carrier  trapping  on  frequency  dispersion  of 
output  resistance  and  transconductance  is  incorporated  in 
Fig.  1.  Following  the  treatment  reported  by  Golio  et  al  [12], 
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Fig.  1.  GaN  MESFET  circuit  model.  Parameters  for  1.5  X  150  /im  GaN 
MESFET  with  Vds  =  25  V  and  Vgs  =  -1  and  -3  V  are  shown  in  the  table. 

the  effects  of  traps  are  incorporated  through  the  parameters 
9m2,Rds2y  and  Cs3.  These  parameters  are  obtained  once  the 
underlying  physical  processes  governing  trap  dynamics  are 
formulated.  Intrinsic  circuit  parameters  Cg^,Cgd,9miiRi,  and 
Rdsi  are  obtained  from  conventional  small-signal  MESFET 
analysis  in  the  absence  of  traps.  Calculated  /-V  characteristics, 
as  shown  in  Fig.  2,  are  based  upon  a  physical  mechanism  for 
current  collapse  and  restoration  proposed  by  the  present  authors 
in  [15].  At  high  drain  bias,  the  electrons  are  captured  by  the 
traps  located  in  the  substrate,  forming  a  second  depletion  layer 
in  the  channel  at  the  channel-substrate  interface  as  shown  in 
the  inset  of  Fig.  2.  This  depletion  layer  gradually  disappears 
with  increasing  drain  bias,  as  the  detrapping  process  is  initiated 
on  the  subsequent  I-V  trace.  The  thicknesses  of  the  depletion 
layer  in  the  channel  under  the  gate  due  to  an  applied  gate  bias 
Vgs  is  h{y)  =  \/{2e/q7ic)[Vbi  +  1^(2/)  +  Vgs]  and  that  at 
the  channel-substrate  interface  due  to  captured  electrons  is 
x„(i/)  =  y/(2e/qnc){l/il  +  nc/NT))[Ubi  +  where 

V\yi  =  {kT/q)\n{nc/ni)  is  the  built-in  potential  of  the 
metal-semiconductor  junction,  17bi  =  (^T/g)  ln(ncJVr/n?)  is 
the  built-in  potential  of  the  channel-substrate  junction,  V{y)  is 
the  potential  along  the  channel  at  y  (i/  =  0  at  the  source)  due  to 
Vbs,  is  the  occupied  trap  concentration  at  the  channel-sub¬ 
strate  interface,  Ui  is  the  intrinsic  carrier  concentration,  e  is 
the  GaN  dielectric  constant,  q  is  the  electronic  charge,  k  is  the 
Boltzmann  constant,  T  is  the  temperature  in  Kelvin,  and  Uc  is 
the  channel  electron  concentration,  ric  equals  the  summation 
of  the  doping  concentration  in  the  channel  Nd  and  the  carriers 
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Fig.  2.  l-V  characteristics  (Vgs  =  and  -3  V)  considering  trapping 
effects.  Calculated  and  measured  results  are  shown  by  solid  lines  and  symbols, 
respectively  [18].  The  inset  shows  the  depletion  regions  in  the  GaN  MESFETs 
when  carriers  are  captured  by  the  substrate  traps. 

recovered  due  to  detrapping  for  a  given  Vds  -  Drain  current  can 
be  expressed  as  Id  =  ~  Hv)  ^n{y)]  where 

fin  is  the  electron  mobility,  d  is  ^e  thickness  of  the  channel, 
W  is  the  width  of  the  device,  and  E{y)  —  dV{y)ldy  is  the 
electric  field  along  the  channel.  Integrating  Id  over  gate  length 
L  (between  y  =  0  and  y  =  I),  the  drain  current  is  expressed  in 
the  following  form: 

Id  =  Ip  {3  {ul  -ul)-2  {ul  -  ug)  -  2[1  +  n, 

•  {VDs)/^T(Vbs)]  {xl  -  x^)}  ,  Vbs  <  VDS.sat  (D 
where 

Ip  ^(^nliXnWd^l&eL 
«d=7(Vbi  +  Vbs-VGs)/Vpi 
«o  =y'(v'bi  -  VGs)/Vpi 
Xd  =  y(i7bi  +  VDsWpi 
xo 

^^Sjsat  ^o) 

Vpi  =qncd‘^/2e 

Upi  =qncd^(l  +  nc/JVr)/2e. 

To  account  for  velocity  saturation  at  the  drain  end  of  the  channel 
at  high  drain  bias,  the  saturation  drain  current  is  expressed  in  the 
following  form: 

Id  =  qNdVsatWd'y{l  -  Vds  >  Vps.sat 


(2) 


84S 

where  7  is  the  saturation  factor,  and  Vsat  is  the  electron  satura¬ 
tion  velocity.  Us  is  obtained  by  equating  drain  current  expres¬ 
sions  of  (1)  and  (2)  at  Vbs  =  ^DS,sat. 

A.  Determination  of  Intrinsic  Model  Parameters 

Intrinsic  circuit  parameters  are  obtained  from  (1)  in  the  ab¬ 
sence  of  trapping  effects  (putting  Uc  =  Nd  and  dropping  the 
last  term  in  the  Vds  <  ^DS,sat  region).  Dividing  the  channel 
into  linear  (length  =  Li)  and  saturation  (length  =  L  -  I/i) 
regions  to  incorporate  channel  length  modulation,  the  intrinsic 
transconductance  Qmi  is  given  by  (3),  shown  at  the  bottom  of  the 
page,  where  Cn  =  (/Xn(l  “  wo))/(7Vsat(l  -  wi))  and  C12  = 
(Li)/(1  -  «i)  -h  {qNdd^^nUi)l{e^Vs^x)-  For  a  given  gate  and 
drain  bias,  Li  and  ui  are  calculated  by  simultaneously  solving 

Li  =  (giVdd2^^)/(6e7Vsat(l-wi))[3(w?-Wo)-'2(wf-ng)] 
andV;,i(?i?”ug)+(2£;crC?^ii)/(7r)sinh[(7r(L-Li))/(2dui)]  = 

Vbs,  which  are  obtained  by  equating  drain  currents  at  L\  and 
solving  Poission’s  equation  in  the  depletion  region.  Intrinsic 
drain  resistance  Rdsi  is  given  by  (4),  shown  at  the  bottom  of 
the  page.  Capacitances  Cgs  and  Cgd  are  given  by 


C^=  —  \ 

(  eqNd'Y  1 

1 

TT 

1 

(5) 

s*  2  ' 

l  Vu) 

2 

1  VLi, 

(6) 

where  Vdg  is  the  drain-to-gate  voltage.  The  second  term 

in  (5) 

and  (6)  represents  the  contributions  from  drain  and  source  side- 
walls  [16]. 

Intrinsic  resistance  is  given  by  jRi  «  (l/^5mi)  where 
T]  «  2[1  -f  (2^Li)/(3(l  -  I3){L  +  Lgd))]^,/5  « 

(l/(  v^l  +  2(Vgs  -  VoFF)tJ'n/Lvsat)  and  Lgd  is  the 
gate-drain  separation  [17].  The  turn-off  voltage  Vqff  is 
estimated  to  be  -6  V. 

B.  Determination  of  Trap-Related  Circuit  Parameters 

The  calculation  of  transconductance  in  the  presence  of  traps 
9m2  requires  an  understanding  of  the  detrapping  processes  at 
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the  channel-substrate  interface.  As  shown  in  Fig.  2,  the  exper¬ 
imental  data  in  the  absence  of  light  show  current  collapse  and 
is  attributed  to  the  presence  of  traps  at  the  channel-substrate  in¬ 
terface.  While  in  the  presence  of  light,  the  trapped  carriers  are 
detrapped  quickly,  and  no  current  collapse  was  observed  [18]. 
It  should  be  mentioned  that  current  collapse  was  not  observed 
in  the  first  measurement  at  Vbs  =  0  V.  However,  current  col¬ 
lapse  was  observed  in  subsequent  measurements  if  the  drain  bias 
was  over  30  V  in  previous  measurements.  As  explained  by  the 
present  authors,  the  observed  characteristics  can  be  explained  by 
identifying  acceptor  traps  with  dual-energy  states  [15].  In  GaN, 
the  impurity  responsible  for  carrier  trapping  that  leads  to  dc  cur¬ 
rent  collapse  has  not  been  identified;  however,  it  is  highly  likely 
that  it  is  due  to  C,  which  may  introduce  multiple  states  within 
the  band  gap,  and  we  refer  to  the  two  states  as  and  X'**'*' 
in  this  analysis.  An  electron  transition  from  the  valence  band  to 
the  neutral  state  leads  to  the  charge  state  X"^*^  (reminis¬ 
cent  of  the  Cr  states  in  GaAs).  This  capture  of  electrons  in  the 
buffer  results  in  the  formation  of  a  negatively  charged  region 
at  the  channel-buffer  interface  consequently  forming  an  addi¬ 
tional  depletion  region  at  the  channel-buffer  interface  causing 
drain  current  to  decrease.  With  increasing  drain  bias,  the  tran¬ 
sition  of  electrons  from  X++  level  to  the  conduction  band  is 
facilitated  as  the  X”*”^  level  moves  closer  to  the  quasi-Fermi 
level.  The  emission  of  electrons  from  the  traps  results  in  a  de¬ 
crease  in  the  width  of  the  depletion  region  at  the  channel-buffer 
interface  with  a  consequent  increase  in  drain  current.  At  high 
drain  bias  (Vds  >  30  V),  carriers  generated  due  to  impact 
ionization  are  trapped,  and  the  trap  states  become  negatively 
charged,  requiring  a  positively  charged  depletion  region  in  the 
channel  at  the  channel-substrate  interface,  and  the  process  con¬ 
tinues.  Using  the  PISCES-II  Device  Modeling  Program  (Ver¬ 
sion  9009-Win32  x  86),  the  simulated  electric  field  exceeded 
1.0  X  10®  V/cm  for  the  structure  reported  by  Binari  et  al  [18] 
for  an  applied  drain  bias  of  30  V.  Kunihiro  et  al  [9]  have  re¬ 
cently  reported  that  in  GaN,  impact  ionization  is  initiated  at 
channel  electric  field  around  0.4  x  10®  V/cm;  therefore,  impact 
ionization  in  the  structures  reported  by  Bianri  et  al  [18]  can  be 
expected.  With  increasing  drain  bias,  beyond  Vds  «  5  V,  the 
quasi-Fermi  level  at  the  drain  end  of  the  channel  comes  closer 
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(3) 


i^dsl  = 


diD 

aVbs 

Vbs 

=  < 


qNdWdfjLn{^-Ud)' 
1 


Vbs  <  ^DS,sat 


qNdVsaaWd 
Vds  >  Vbs, sat 


1 2wi  Vpi  -h  sinh 


7r(L-Li) 
2dui 


—  Ect 


^12  +  - 


- — cosh 

Ui  ) 


7r(L-Li) 


2dui 


]} 


(4) 


ISLAM  et  al.:  PHYSICS-BASED  FREQUENCY  DISPERSION  MODEL  OF  GaN  MESFETs 


849 


to  the  trap  level,  thereby  increasing  the  cross  section  of  elec¬ 
tron  transition  from  the  trap  level  to  the  conduction  band.  The 
initiation  of  electron  detrapping  reduces  the  width  of  the  deple¬ 
tion  region  that  results  in  an  increase  in  drain  current.  Beyond 
Vbs  «  30  V,  impact  ionization  is  initiated;  the  electrons  are  cap¬ 
tured  by  the  acceptor  traps,  and  the  process  continues.  The  drain 
voltage  at  the  onset  of  detrapping  process  (~5  V  at  Vbs  =  -1 
V,  from  Fig.  2)  is  denoted  as  which  decreases  with 

decreasing  Vbs  as  evident  from  Fig.  2.  VjQNSET  obtained  by 
solving  ( Vbi + Vgj^^^'^-yQs)^  1  /  ( V 1  +  Nd/Nto)  ( t/bi  -h 

^^NSET^i/2  _  which  is  obtained  by  applying 

the  pinch-off  condition  to  (5/D)/(^Vbs)|  Vbs  =  facilitate 
calculation,  the  computed  Vbs  dependence  of  is  ap¬ 
proximated  by  the  quadratic  form  =  ao  +  oiVbs  + 

^2^03  where  ao  =  6.24  V  and  ai  =  1.29  and  02  =  0.03  V”^. 

Theoretical  calculations  based  upon  the  discussed  model  are 
compared  with  experimental  data,  in  Fig.  2,  to  show  good  agree¬ 
ment.  Calculations  are  carried  out  for  a  1.5  ^m  x  150  /xm  GaN 
MESFET  with  2000- A  n-GaN  channel  grown  over  a  3  /xm  semi- 
insulating  (SI)  GaN  substrate.  The  channel  doping  concentra¬ 
tion  JVd  is  2  X  10^*^  cm”^.  The  measured  low-field  mobility 
is  410  cm^  •  •  s“^  [7],  [18].  The  calculation  is  based  upon 

the  assumption  that  the  drain  bias  dependence  of  trapped  earner 
concentration,  in  the  drain  bias  range  between  and  the 

drain  bias  at  the  onset  of  impact  ionization  (^^30  V  at  Vqs  =  —  1 
V),  is  Nt  =  Nto  exp[-a(FDS  -  At  the  gate  end  of 

the  drain  the  current  due  to  the  detrapping  of  carriers  is  mod- 
eledas/Detrap  =  9«satiVtoW"d[l-exp{-a(Vbs-V;?sNSET)}] 
where  Ugat  is  the  carrier  saturation  velocity  in  the  channel.  The 
fit  with  the  measured  I-V  characteristics  as  shown  in  Fig.  2,  is 
obtained  with  Nto  =  1.5  x  10^®  cm”"^  and  a  =  0.112  V“^. 
The  transconductance  representing  the  effect  of  detrapping  is 
obtained  by  differentiating  the  detrapping  current  with  respect 
to  I^dg  as  shown  below 


5m2  =  qVsB^t^dNtod 


ai  +  2a2Vbs 
1  -f  ai  -f-  2a2Vbs 


xexp[-a{VDs-V^s^^^)]-  O) 


By  defining  W(,(rds)  as  u»<,(rds)  =  !/*<«  where  Wo(rds)  is  the 
frequency  at  which  overall  output  resistance  has  decreased  to 
the  average  of  its  low-frequency  and  the  high-frequency  values 
Eds2  is  given  by. 


R<Is2  «  ^  B  V 

G^8s(1  +  OmlRdsl) 

where  td  is  the  detrapping  time  constant.  Css  is  comparable 
to  Cgs  [12]  and  the  other  circuit  parameters  are  obtained  from 
reported  experimental  data  and  are  as  follows:  =  6  0,  = 

70  U,Rs  =  90f2,and(7ds  =  0.040  pF  [18].  However,  due  to  the 
unavailability  of  data  for  GaN-based  devices,  the  inductances 
associated  with  connecting  pads  are  assumed  to  be  the  same 
as  those  of  GaAs  FETs  and  are  as  follows:  Lg  =  0.055  nH, 
Ld  =  0.307  nH,  and  Ls  =  0.027  nH  [12]. 

The  determination  of  the  detrapping  time  depends  upon  the 
nature  and  location  of  the  traps.  An  estimation  of  td  based  upon 
the  drain  lag  measurement  reported  by  Binari  et  al.  [8],  for 


Fig.  3.  Extracted  detrapping  current  (symbols)  from  the  drain  lag  measure¬ 
ment  reported  in  [8].  Solid  lines  show  the  fits  using  single  and  multiple 
detrapping  time  constants  [23].  Output  resistance  dispersion  frequency 
(/o(rds))  as  a  function  of  temperature  is  shown  in  the  inset  for  different  trap 
levels  involved  in  the  detrapping  process. 


AlGaN-GaN  HEMTs,  is  shown  in  Fig.  3.  Assuming  exponen¬ 
tial  decay  of  the  detrapping  current  the  detrapping  time  con¬ 
stant  is  estimated  to  be  53.2  s  and  is  in  agreement  with  the 
value  mentioned  in  [8].  The  detrapping  time  is  expressed  as 
[11]:  fd  «  l/(tTnt^thiVcexph(£;c-^t)/fcT]},  where  an  is  the 
electron  capture  cross  section,  uth  =  \/3kT/m^  is  the  thermal 
velocity,  m*  is  the  electron  effective  mass  in  GaN,  Ec  and  Et 
are  the  conduction  and  trap  energy  levels,  respectively.  With 
m*  =  0.2mo,  where  rrio  is  free  electron  mass,  iVc  =  1  x  10^^ 
cm“^,  and  tT„  =  cm“^,J5t  is  found  to  be  located  at 

0.78  eV  below  Ec  at  room  temperature.  A  better  fit  to  the  detrap¬ 
ping  current,  at  300  K,  is  obtained  by  using  multiple  time  con¬ 
stants  as  shown  in  Fig.  3.  The  time  constants  of  58.42  and  1.55  s 
correspond  to  trap  states  located  at  0.79  and  0.69  eV  below  Ec, 
respectively.  As  observed,  these  trap  levels  at  room  temperature 
correspond  to  dispersion  frequencies  of  less  than  1  Hz.  How¬ 
ever,  with  increasing  temperature  the  dispersion  frequency,  due 
to  the  traps  located  at  0.79  and  0.69  eV,  increases  as  shown  in 
the  inset  and  can  be  of  the  order  of  MHz  at  600  K.  Moreover, 
a  higher  dispersion  frequency  is  possible  due  to  the  presence  of 
other  shallow  traps.  Although  the  estimated  trap  levels  are  based 
upon  data  reported  for  an  AlGaN-GaN  HEMT  the  calculated 
trap  levels  are  consistent  with  the  reported  experimental  values 
for  GaN.  Reported  trap  levels  measured  from  the  conduction 
band  minimum  are  0.59  and  0.66  eV  [19],  0.598  and  0.67  eV 
[20],  0.59  and  0.88  eV  [21],  and  0.41-0.50  and  0.62-0.67  eV 
[22]. 

T-parameters,  as  listed  in  the  Appendix,  are  calculated  by 
placing  short  circuits  at  the  output  and  input  terminals  of  the 
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Fig.  4.  Transconductance  as  function  frequency  with  Vbs  as  parameter 
(Vbs  =  —1  V  and  detrapping  time,  f  j  =0.5  ms)  [23], 

circuit  shown  in  Fig.  1. 121  and  Y22  are  simplified  at  low  and 
high  frequencies  to  formulate  the  overall  transconductance  and 
output  resistance. 

m.  Results  AND  Discussion 

Fig.  4  shows  the  variation  of  the  overall  device  transcon¬ 
ductance  as  a  function  frequency  with  drain  voltage  as 
parameter.  The  gate  voltage  is  —1  V  and  the  detrapping 
time  is  assumed  to  be  0,5  ms.  For  frequencies  lower  than 
fo(gm)  «  l/2wRds20ss,  wWch  is  the  dispersion  frequency 
for  transconductance  [12],  the  overall  transconductance  equals 

=  (5mi/l+5mii?s+(i^d+i^3/i^dsi))andforfrequencies 

greater  than  fo(gm)  the  overall  transconductance  equals  « 
(dml  ~  ffmlRs  4" gm2Rd  (Rd  ~i~  Rs) / (-^dsl  ||  ^ds2))* 

The  above  expressions  are  obtained  by  simplifying  Y21.  The 
overall  transconductance  increases  with  increasing  Vos  until 
the  saturation  voltage  is  reached.  The  increase  in  the  overall 
transconductance  at  low  frequency  is  due  to  the  increase  in  gmi 
with  increasing  Vbs  and  the  decrease  at  higher  frequencies,  for 
a  given  Vbs»  is  due  to  the  contribution  of  gm2  as  stated  earlier. 
At  elevated  drain  bias,  most  of  the  trapped  electrons  have 
already  been  released  and  the  detrapping  current  is  very  small 
that  results  in  a  lower  transconductance  dispersion,  which  is 
also  evident  from  (7). 

Fig.  5  shows  the  frequency  dispersion  of  output  resis¬ 
tance  with  drain  bias  as  parameter  at  Vgs  =  -1  V.  For 
finite  Rs  and  Rd  the  overall  output  resistance  R^'^  = 


Fig.  5.  Output  resistance  as  function  frequency  with  Vbs  as  parameter 
V  and  detrapping  time,  f<£  =0.5  ms).  Inset  shows  the  variation  of 
Rdsi,Rds2*  and  gm2  as  a  function  of  Vbs  for  Vbs  =  —1  V  [23]. 

Rs  Rd  RdsiO-  +  9miRs)  for  frequencies  less  than  the 
output  resistance  dispersion  frequency  /o(rds)  ^  l/27ctd^  For 
frequencies  greater  than  /o(rds)»  the  overall  output  resistance 
R^^  «  (i^dsl||i^ds2||l/^m2)(l  +  9mlRs  +  9m2Rd  +  {Rd  + 
Rs)/{Rdsi  Il-Rds2))-  The  above  expressions  are  obtained  by 
simplifying  Y22-  At  low  Vbs,  the  output  resistance  does  not 
show  any  dispersion  and  equals  For  frequencies  lower 
than  /o(rds)  «  l/27rfd,  only  Rdsi  contributes  to  the  overall 
resistance  as  the  branch  containing  Css  behaves  like  an  open 
circuit  due  to  its  high  capacitive  reactance.  Moreover,  with 
increasing  Vos  the  slope  of  the  current- voltage  characteristics 
decrease  resulting  in  the  observed  increase  in  R^^-  It  is  seen 
from  (8)  that  Rds2  is  always  much  larger  than  Rdsi  and  l/gm2 
(see  inset  of  Fig.  5)  and  for  fi*equencies  greater  than  /o(rds)  the 
overall  output  resistance  is  mainly  determined  by  the  parallel 
combination  of  i^dsi  and  l/pm2*  With  increasing  drain  bias 
Rdsi  increases  and  due  to  decreasing  detrapping  current  gm2 
decreases  (see  inset  of  Fig.  5).  Therefore,  dispersion  of  the 
output  resistance  is  maximized  at  a  certain  drain  bias  beyond 
which  dispersion  decreases  and  the  overall  output  resistance 
approaches  R^^  and  is  attributed  to  the  diminishing  values 
of  5m2-  A  similar  behavior  was  observed  in  output  resistance 
dispersion  in  GaAs  FETs  [13],  [14]. 

Fig.  6  shows  the  dispersion  in  transconductance  as  a  function 
ofgatebias.  The  drain  bias  is  10  V  and  detrapping  time  is  0.5  ms. 
The  overall  low-frequency  transconductance  is  governed 
by  the  behavior  of  gmi  which  increases  with  increasing  gate 
bias.  The  saturation  of  gmi  causes  the  observed  rate  of  increase 
of  g^  with  Vgs  to  decrease  at  higher  gate  bias.  At  higher  fre¬ 
quencies,  the  overall  transconductance  dispersion  g]^  -  g^  is 
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Frequency  (Hz) 

Fig.  6.  Transconductance  as  function  ftequency  with  Vbs  as  parameter 
(Vbs  =  10  V  and  detrapping  time,  td  =  0.5  ms). 


With  decreasing  gate  bias,  Rdsi  increases  and  the  parallel  com¬ 
bination  of  ildsi,  l/9m2  and  Rds2  is  dictated  by  l/gm2  resulting 
in  a  higher  dispersion. 

IV.  Conclusion 

Frequency  dispersion  of  transconductance  and  output  resis¬ 
tance  in  GaN  MESFETs  is  reported  by  considering  the  effect  of 
traps.  Detrapping  time  is  estimated  from  drain-lag  measurement 
data  reported  for  AlGaN-GaN  HEMT.  The  trap  levels  located 
at  0.69  and  0.79  eV  are  comparable  to  the  experimental  data 
reported  for  bulk  GaN.  The  dispersion  frequencies  estimated 
from  the  trap  levels  are  in  the  range  of  MHz  at  elevated  tempera¬ 
tures,  where  a  typical  GaN  power  device  may  operate,  although 
at  room  temperature  it  may  be  few  Hz.  In  the  linear  region  of  op¬ 
eration,  dispersion  in  transconductance  is  significant,  whereas, 
dispersion  in  output  resistance  is  crucial  when  the  device  oper¬ 
ates  near  cutoff. 

Appendix 

y-parameters  of  the  circuit  shown  in  Fig.  1,  obtained  by 
placing  short  circuits  at  the  output  and  input  terminals  are  as 
shown  at  the  top  of  the  next  page,  where 


I.OE‘1'00  1.0E+02  1.0E+04  1.0E+06 

Frequency  (Hz) 

Fig.  7.  Output  resistance  as  function  frequency  with  Vbs  as  parameter 
(Vbs  =  10  V  and  detrapping  time,  td  =  0.5  ms). 


observed  to  decrease  with  increasing  Vqs-  This  decrease  is  at¬ 
tributed  to  the  effect  of  different  resistive  components  in  the  cir¬ 
cuit  as  evident  from  the  expressions  of  and  g};!^ .  It  should  be 

noted  that  with  increasing  Vgs  although  -  g^  decreases, 
gm2  increases  slightly  due  to  an  increase  in  Fig.  7 

shows  the  dispersion  of  output  resistance  with  varying  gate  bias. 


^1  —  TJni 

^2  “  ^sl 

=^Yg  y^s  +  Yml  +  Tdsl  +  Yss 

AI4  ”  Ygg 

hi  —  ^^d  Yjnl  Yjfi2 

h2  =  Ygd  +  Idsl  +  Ym2  +  yds2  +  Yd 
hz  =  —y^ml  —  Yds! 
h^  —  yds2 

Ai  =  (ygs  +  n  +  ygd)/rgs 

YssiY^s  +  ~  ym2ygs 

ygs(n.+ids2) 

Bi  =  —Ygi/Ygs 

ygB(yds2+ym2)-ygdy.. 

^  Fgs(n.  +  l^ds2) 

Cl  =  -Y,/Yg, 

r  nn, 

ygs(i;«  +  i'ds2) 
yml  =  5Tnl/(l  +i<^C'gsiii) 
im2  —  9m2 
Yg  =  l/{Rg  +  jU)Lg) 
is  =  1/(^S  + 

Yd  =  l/{Rd  +  ji^Rd) 

igs  “  J^Cgs/(l  "i  j^^gsRi) 
igd  =  iwCgd 

Ys8  “ 

idsl  —  (1  +  J^CdsJ?dsl)/'^sl 
ids2  =  1/R<is2- 
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Vn  =  y, 
Yi2  = 


1  + 


{k^Ci  +  k4C2){h2  +  hzBi  +  /i4-^2)  ^  (^3^1  +  h^C2){k2  +  ^3^1  ~t~  ^4-^2)  1 

(All  4-  k^Ai  +  k4A2){h2  +  h^Bi  +  /i4^2)  “  (A^2  +  k^Bi  +  k4B2){hi  +  h^Ai  +  /i4^2)  J 


{hi  +  h^Ai  +  /i4-^2) 

{k2  +  k^Bi  4-  A:4J92)(^i  H"  ^3^1  H"  ^4-^2)  ”  (A^i  “I"  k^Ai  4*  k4A2){h2  +  ^3-^1  4-  ^4-^2)] 


(/i2  4-  hsBi  4-  h4B2){ki  4-  k^Ai  4-  ^4^2) 


1^21  = 


-Yd 


—  {h^Ci  4-  ^4^2)  4"  {hi  4*  h^Ai  4-  ^14^2) 


{h2  4-  h^Bi  4-  ^4-^2) 

(A:3C'i  4-  k4C2){h2  +  h^Bi  4-  /t4-^2)  "  (^3^1  4-  ^4^2)(fe2  4-  ^3^1  4-  ^4^2) 


1^22  =  Yd 


{k2  4-  ^3^1  4-  A:4^2)(^2  4-  hzBi  4-  /i4-S2)  -  (A:2  4-  kzBi  4-  k4B2){hi  4-  hzAi  4-  /i4^2) 

^d(A:i  4*  ^3-^1  4“  ^4^2) _ 


1  + 


(/c2  4-  kzBi  4-  k4B2){hi  4-  hzAi  4-  /i4^2)  —  (A^i  4-  A:3i4i  4-  k4A2){h2  4-  hzBi  4-  ^14^2) 
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